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A field theory of the Anderson transition in two dimensional disordered systems with spin-orbit 
interactions and time-reversal symmetry is developed, in which the proliferation of vortex-like topo- 
logical defects is essential for localization. The sign of vortex fugacity determines the Z2 topological 
class of the localized phase. There are two distinct, but equivalent transitions between the metallic 
phase and the two insulating phases. The critical conductivity and correlation length exponent of 
these transitions are computed in a N = 1 — e expansion in the number of replicas, where for small 
e the critical points are perturbatively connected to the Kosterlitz Thouless critical point. Delo- 
calized states, which arise at the surface of weak topological insulators and topological crystalline 
insulators, occur because vortex proliferation is forbidden due to the presence of symmetries that 
are violated by disorder, but are restored by disorder averaging. 
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Topology can have a profound impact on Anderson lo- 
calization in disordered electronic systems. This was first 
understood in the integer quantum Hall effect [IJ [2] , where 
the two dimensional (2D) bulk states at the plateau tran- 
sition are extended, even in the presence of strong dis- 
order. In the field theory of localization, this dereal- 
ization is associated with the presence of a topological 
term in the nonlinear sigma model (NLctM)[3]- At long 
distances, generically renormalizes to integer multiples 
of 27T reflecting the topological quantization of the Hall 
conductivity, while = -k is a critical point, at which 
bulk states are extended. More recently, this paradigm 
for derealization has been extended to other symmetry 
classes and dimensions [U [5]. 

It is now known that there are two topologically dis- 
tinct classes of 2D insulators in the presence of time rever- 
sal (TR) symmetry [6 --8J. A field theory of localization in 
the symplectic class should incorporate and distinguish 
the two insulating states. Such a theory must go be- 
yond the single-parameter scaling based on conductivity 
alone [9HI2]. but it has been unclear what additional pa- 
rameter plays the role of 9. 

A related difficulty is revealed by recent studies of sur- 
face states of 3D weak topological insulators (WTI) [TBT 
[T3] and of topological crystalline insulators (TCI)[Tt)l 
117] . These surfaces have an even number of 2D Dirac 
fermions, so the Zi topological term that guarantees de- 
localization of a 3D strong topological insulator (STI) 
surface [3 HH] is absent. Nonetheless, general arguments, 
as well as numerics [13], suggest that contrary to the con- 
ventional picture, the surface remains delocalized even 
with strong disorder, due to the existence of discrete sym- 
metries that are violated by disorder, but remain unbro- 
ken on average. This poses the question of how average 
symmetries fit into the field theory of localization, and 



how those symmetries prevent localization. 

In this paper we answer those questions by examin- 
ing the crucial and unexplored role played by topologi- 
cal defects in the NLcrM in the 2D symplectic class|19j. 
We show that the proliferation of pointlike vortices is 
essential for localization, and that the sign in the parti- 
tion function associated with those vortices distinguishes 
a topological insulator (TI) from a trivial insulator in 
2D. The ensemble average symmetry forces the system 
to be on a line where the fugacity of vortices vanishes 
and hence dictates derealization. In addition to explain- 
ing the origin of the metallic state, this analysis provides 
new insight into the 2D symplectic metal-insulator tran- 
sition. We show that there are two distinct but equivalent 
fixed points that describe transitions to insulator and TI 
states. Moreover, by treating the number of replicas, N 
as a continuous variable, we show that for N = 1 — e the 
Anderson transition fixed points are perturbatively con- 
nected to the Kosterlitz- Thouless (KT) transition fixed 
point [2D] for e — > 0. This allows us to compute the critical 
conductivity and correlation length exponent perturba- 
tively in an e expansion. 

Before describing the symplectic class, we briefly dis- 
cuss a simpler version of derealization via average sym- 
metry in the unitary class. It is well known that the 
surface of a TR invariant 3D STI is delocalized[THl |2~T] . 
but TR violating perturbations can lead to localization. 
This localized state is in a sense "half of a quantum Hall 
state, and has a xy — ±e 2 /2h. Importantly, the time re- 
verse of this state, with a xy — =Fe 2 /2/i, is topologically 
distinct. If impurities have random local moments, TR 
symmetry remains unbroken on average. This places the 
system precisely at the transition between the two local- 
ized states. This can be modeled by considering an en- 
semble in which each member violates TR, but the whole 
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ensemble is TR invariant. Such a system is described by a 
NLerM in the unitary class, which in 2D allows a Pruisken 
6>-term. Since TR, when applied to the ensemble, takes 9 
to —6, average TR symmetry constrains 9 to be or tt. 
The surface of a STI corresponds to 9 = tt, so the system 
flows to the delocalized critical point associated with the 
plateau transition. If the average symmetry is broken 
by an applied magnetic field, then the system flows to a 
localized phase with a xy = ±e 2 /2h. It would be inter- 
esting to probe this experimentally, and to compare the 
conductivity and critical behavior as a function of applied 
magnetic field with integer quantum Hall transitions. 

The surface of a WTI or a TCI also possesses an aver- 
age symmetry. For the WTI, the symmetry is translation 
by one lattice constant in the direction of the layering of 
the WTI, while for the TCI studied in Ref. H3 it is a 
mirror symmetry. By breaking this symmetry, it is possi- 
ble to gap the surface, leading to localization. However, 
applying the discrete symmetry to the gapped state leads 
to a topologically distinct localized state. Interfaces be- 
tween the two localized states are associated with a heli- 
cal edge mode. If the discrete symmetry is respected on 
the average, then the system is precisely on the bound- 
ary between the two localized states — analogous to the 
transition between the 2D TI and trivial insulator. Un- 
like the unitary case, this transition need not occur at a 
point, but can involve an intermediate metallic phase. It 
is clear, however, that even for strong disorder, the sys- 
tem can not be localized at this point because it is im- 
possible to change the topological class without extended 
states crossing the Fermi energy. 

To develop a field theory for this derealization, we 
use the fermionic replica theory introduced by Efetov 
et al.jTTj. Our analysis closely parallels that of Ryu, et 
al.[ISJ. We consider a system with average Hamiltonian 
"Ho j along with gaussian correlated TR invariant disor- 
der. Using the replica trick to integrate out the disorder, 
the disorder averaged product of retarded and advanced 
Green's functions can be generated from the partition 
function Z = J D[^p,Tp]e~ s , with 

S = J d 2 r[tp a ((7io-E)5 ab + ir]A ab )ij; b - ^{i) a '4 ! b){^btpa)] 

(1) 

Here a = 1, ...,2N is an index for N retarded and N 
advanced replicas, and A = ljy © (— 1)at, where 1 at is a 
N x N identity matrix. ip a is a Grassman field, which in- 
cludes (suppressed) spin, position and possibly orbital in- 
dices. ip a = ipai<J v , where a acts on the spin indices. TR 
symmetry requires a v V.ocr v = %o, so that ia v (H — E) is 
a skew symmetric matrix. For rj = 0, |T]) is invariant un- 
der 0(2N) rotations among the replicas, which is broken 
down to O(N) x O(N) by r). 

A theory of the Nambu Goldstone modes associated 
with this symmetry breaking is formulated by Hubbard 
Stratonovich decoupling the four fermion interaction, and 



performing a saddle point expansion about the broken 
symmetry state. After freezing the massive modes, the 
saddle point is characterized by a 2N x 2N matrix field 
Q = T KO, with O £ 0(2N). Distinct values of Q 
belong to the coset G/H = 0(2N)/0(N) x O(N) and 
satisfy Q = Q T ', Q 2 = 1. A theory for the long wave- 
length fluctuations in Q ab is obtained by integrating ip a 
in the background of a spatially varying Q ab . This gives 
Z cS = J D[Q]e~ s ^^ with 

e -S off [Q] _ J £)^p ; ^ e -J d 2 r[^ a [{H -E)8 ab +i^Q ab \^ b \ 

(2) 

Here A is a parameter characterizing the bare scattering 
time that is determined self consistently at the saddle 
point. Expanding in gradients gives the NLcM, 

S° eS [Q] = 3^/ d 2 rTr[(VQ) 2 ], (3) 

where the coupling constant t characterizes the disorder 
strength and is related at lowest order to the resistivity, 
a = (2wt)~ 1 e 2 /h. The renormalization of t at long wave- 
lengths is described by the perturbative renormalization 
group (RG) equation[TTJ HlllUll] 

dt/d£ = (3(t), (3(t) = 2(N - l)t 2 + ... (4) 

In the replica limit, N — > 0, the weak coupling fixed point 
t = is stable, indicating the stability of the symplectic 
metal phase, characterized by weak antilocalization. 

Eq.([3| is not the whole story because topologically 
non trivial configurations of Q can have important non- 
perturbative effects. There are two types of topologi- 
cal configurations associated with the first and second 
homotopy group of G/H. Ryu et aZ.[18] argued that 
tt 2 (G/H) = Z 2 allows for a topological term with 9 = n 
associated with the Z 2 global anomaly. Such a term ap- 
plies to the surface states of 3D STIs, but is absent in 
2D systems as well as the WTI or TCI surface because 
of fermion doubling. 

For our problem, the crucial topological objects are 
point-like defects in 2D similar to vortices which are al- 
lowed by the nontrivial ir\{G/H) — Z 2 - We now argue 
that these defects are necessary for localization, and that 
their contribution to Z e g encodes the distinction between 
a trivial and topological insulator. The role of vortices 
can be understood by considering an inhomogeneous 2D 
system in which a TI in region S with boundary C is 
surrounded by a trivial insulator, as shown in Fig. 1(a). 
Consider the effect of integrating out tp a in ^ for weak 
disorder in the presence of a vortex configuration Q ab (r) . 
Since the interior of S has a finite band gap, dominant 
contribution to the action comes from the helical edge 
states at the boundary C. On C, Q a b{ Y € C) is a non- 
singular and non-contractible configuration correspond- 
ing to the nontrivial element of tc±(G/H). Repeating the 
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FIG. 1: (a) A 2D TI (m < 0) in region S with boundary 
C is surrounded by trivial insulator (m > 0). The sign of 
the partition function depends on the parity of the number 
of vortices in S. (b) For m — the eigenvalues of 'H eff in 
( 10 1 exhibit a linear zero crossing, which leads to a vanishing 
fugacity for vortices. 



analysis of Ryu et al. [18] in ID, one can explicitly ver- 
ify that the derealization of the helical edge states is 
associated with a topological term in the ID NLtrM[S], 
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n(C) 
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where n(C) = 0, 1 is the Z% homotopy class of Q on C. 
Importantly, since Q is denned in all space (except at the 
cores of vortices), n(C) may be viewed as the number of 
vortices in S mod 2. This leads to a bulk characterization 
of the TI based on the 2D NLcM: in the TI the fugacity 
v of Z2 vortices is negative. In the trivial insulator, the 
topological term is absent, and v is positive. The posi- 
tivity of v for a trivial insulator can also be understood 
in the limit of vanishing spin-orbit coupling in Hq: the 
fermion determinant from integrating out ip a is squared 
due to spin degeneracy. 

At the transition between the trivial and TI, v must 
pass through zero. This suggests that at the WTI sur- 
face, v = 0. We will now proceed to demonstrate this 
explicitly. To model the WTI surface, consider a sys- 
tem of helical edge modes H = va x k x that are stacked 
in the y direction with separation a. Coupling between 
neighboring layers gaps the surface, except at two Dirac 
points at (k x ,k y ) = (0,0) and (0, 7r/a). If we index the 
two Dirac points by t z = ±1, then the surface states are 
described by the Dirac Hamiltonian, 



H = v(a x k x + a y r z ) + ma y T y 



(6) 



The discrete translation symmetry of the WTI corre- 
sponds to translating by ay, described by exp(ipj,a) = t z . 
We have included a uniform mass term mg^llij, that 
describes a translation symmetry breaking dimerization 
of the layers. This is the only mass term that respects 
TR symmetry. The two topologically distinct dimeriza- 
tion patterns are distinguished by the sign of to. (In- 
cluding a large k regularization a y T v k 2 , this model also 
describes the 2D transition between a trivial (to > 0) and 
topological (to < 0) insulator p5].) 

In general, the NLctM ^ should include a sum over 
vortex configurations in Q. The fugacity of the vortices is 
determined by comparing ([2]) in the presence and absence 



of vortices. To evaluate ([2]) for a vortex consider the 
simplest vortex configuration involving a single retarded 
and advanced pair of replicas. This can be expressed in 
terms of a one parameter family of Q's of the form 
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A Z2 vortex is then a configuration where 9 winds by an 
odd multiple of 27r. 

The Grassman integral in ^ defines a Pfaffian, so that 
the vortex fugacity may be written 



= Pf[igVD{Q)] 
V ~ H[iovD(Q )]' 

Where Q is a vortex configuration, and Qo 
space of the two nontrivial replicas we have 



(8) 



A. In the 



D(Q) = {H - E) +iA(/i 2 cose» + M K sin6»). (9) 

Here fj, z is a Pauli matrix in the space of the two non- 
trivial replicas. To evaluate the Pfaffian, we use a trick 
similar to that used by Ryu et al.l8\, and compute 
(Pi[ia y D}) 2 = det[icr^L>] = det\p y D]. This is useful be- 
cause fi y D = T-L is a Hermitian operator given by 
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H y {U -E) + A(fi x cos 9 - n z sva6), (10) 



so the determinant is the product of its real eigenvalues. 
The TR symmetry of the original Hq becomes a particle- 
hole symmetry, {H cff ,S} = 0, with 3 = [i y a y K. More- 
over, when to = 0, T~L cS decouples into two independent 
Hamiltonians for r z = ±1. Each is identical to a topolog- 
ical superconductor in class D, with 9 playing the role of 
the superconducting phase. It follows that there are two 
zero modes indexed by t z = ±1 bound to the core of a 
vortex with odd-vorticity. For to 7^ 0, the zero modes are 
coupled and split, as shown in Fig. lb. We conclude that 
det{[i y D] has a second order zero at to = 0, so Pi[ia v D] 
has a first order zero, which involves a sign change as a 
function of m. This shows that sign of the fugacity v 
of the vortices depends on to and vanishes for m = 0. 
Thus isolated Z2 vortices are forbidden at the surface of 
a WTI. More generally, it is possible to have multiple 
vortices, since in that case the zero modes will split even 
for m — 0, leading to a nonzero Pfaffian. However, since 
the splitting vanishes exponentially in the separation, the 
vortices will be confined by a linear potential. 

It is thus clear that the vortex fugacity v is a crucial 
variable that must be included in the NLctM. The topo- 
logical and trivial insulators are distinguished by the sign 
of v. Proliferation of vortices leads to a disordered state 
in the corresponding topological class. For v — 0, quali- 
tatively different behavior is expected to occur reflecting 
the derealization of the WTI or TCI surface states. In 
this case, the target space of the NLcM effectively lifts 
to its double cover, G/H = SO(2N)/SO{N) x SO(N), 



for which n^G/H) = 0. Since G/H and G/H have iden- 
tical local structure, their perturbative (3 functions will 
be identical. It is useful to consider this behavior as a 
function of the replica number, N . 

For N > 1, (3(t) > 0, and the weak coupling fixed point 
is unstable, leading to a disordered phase even in the ab- 
sence of vortices. This phase is "less disordered" than the 
disordered phase with v ^ 0, though. The confinement 
of Z-2 vortices leads to a topological order similar to a Z2 
spin liquid |26j . This can be clearly seen by placing the 
system in a torus: there are four topologically discon- 
nected sectors corresponding to the homotopy classes of 
Q G G/H along the two large loops. When v is turned 
on in this disordered phase, the Z2 vortices immediately 
condense. The v — line thus describes a first order 
transition between the v > and v < phases. 

The behavior for N — > is expected to be qualitatively 
different. In this case the weak coupling fixed point is 
stable, and we know from the arguments presented above 
that even at strong coupling there must be derealization 
for v — 0. It is useful to consider the critical value N = 1 
that separates these behaviors. The theory for TV = 1 is 
simply the XY model, and Q is fully parameterized by 9 
in The action ^ becomes 
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Since the target space, S 1 , is flat, f3(t) = to all orders, 
but vortices modify the behavior. For small t, 2it vor- 
tices in 9 are bound, and the system flows to a fixed line 
parameterized by t. For t > t* = 1/16 vortices unbind at 
a KT transitionf2Uj to a disordered phase. 

We now consider the behavior for N < 1, treating N 
as a continuous variable. Since Z2 vortices are present 
for all N, it is reasonable to examine their effects as a 
function on N. We find that the theory can be controlled 
for N — 1 — e, with e«l. To lowest order in e, the KT 
flow equations to lowest order in v are modified by the 
nonzero (but small) f3{t) = (N - l)P(t), 



dt/dl 
dv/d£ 



-ej3(t) + v 2 
(2-(8i)->. 



(12) 



To this order, we are free to set the coefficient of v 2 to 
one by rescaling v. The RG flows are shown in Fig. 2. 
There are two fixed points at 



t* = l/W, v* = ±[ep{t*)] 1 / 2 



(13) 



For small e, these fixed points are within perturbative 
range of the KT fixed point. These fixed points describe 
a transition between the ordered and disordered phases 
of the 0(2N)/0(N) x O(N) NLctM for N < 1. For 
iV^Owe identify these fixed points with the Anderson 
transition between the symplectic metal and the localized 
trivial/topological insulator. These two metal-insulator 
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FIG. 2: (a) RG flow diagram based on (12 1. The stable fixed 
point at (t, v) = (0,0) is the symplectic metal (SM). The two 
unstable fixed points at (t* ,±u*) approach the KT transition 
for e = 1 - iV < 1 and for TV — > are identified with the 
Anderson transition, (c) is a flow diagram which includes an 
third fixed point at (t m , 0), along with a fixed point at (t a , 0) 
which describes a direct transition between TI and I. (b) and 
(d) are phase diagrams corresponding to (a) and (c). 



transitions have identical bulk critical behaviors, which 
is expected from the fact that the total number of Z 2 
vortices in a closed system is always even. 



By expanding (12 1 about the fixed point, we can iden- 



tify the critical conductivity and the correlation length 
exponent associated with the symplectic Anderson tran- 
sition. To lowest order in e we find a* = (2irt*)~ 1 e 2 jh — 
{8/n)e 2 /h, and v = 2t* / {^{t*)) 1 / 2 . While (3(t*) is not 
known exactly, f3(t) has been computed perturbatively 
up to order t 5 [23]. The small value of t* is well within 
the range of this perturbation theory. The second order 
term gives only 6% correction and the higher terms are 
even smaller. Using the first term from Eq. [4] we find 
v = (2/e) 1 / 2 . Extrapolating to e = 1 gives 



a* - 2.5e 2 /h, 



1.4. 



(14) 



These values are rather different from numerical es- 
timates of critical exponents in previous model studies, 
which give a* - lAe 2 /h and v - 2.7H7HH]. We sug- 
gest two possible origins of the discrepancy, depending 
on the behavior of the N = NLerM at strong coupling, 
which cannot be accessed in the present analysis. One 
possibility is that for N — > 0, f3(t) < for all t along the 
line v = 0. The corresponding RG flow and phase dia- 
grams are shown in Fig.2a-b. In this case, the symplectic 
metal-insulator transition is governed by the fixed point 
(t* , v*). The discrepancy in exponents is then most likely 
due to the slow convergence of the e expansion similar to 
the d = 2 + e expansion for the 3D Anderson transition. 
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A more interesting possibility is that for N — > 0, j3(t) 
changes sign at a critical point t m on the line v — 0. In 
fact, t m is present for N = 1 — e. For N = 1, double 
vortices are allowed, and will in general have non zero 
fugacity. The theory with both single and double vortices 
can be analyzed using a dual sine-gordon theory, 

S = J d 2 r-(V(^) 2 +v costp + v 2 cos2ip (15) 

where vi is the fugacity for double vortices. When v = 0, 
V2 becomes relevant at t m — 1/4. When w 2 flows to 
strong coupling, v = describes a first order transition 
similar to the case when N > 1. However, it is un- 
likely that this first order transition persists to N = 0, 
which is a theory of disordered non-interacting electrons. 
Instead, the most likely scenario is that there is a con- 
tinuous direct transition between trivial and topological 
insulators, controlled by a strong coupling fixed point t s , 
as indicated in Fig. 2c-d. If this is the case, the flow 
trajectory near v = line may be strongly influenced 
by the two nearby fixed points t m and t s . Nonetheless, 
the true scaling behavior of metal-insulator transition in 
the symplectic class is ultimately controlled by the fixed 
points (t*,v*). In this regard, we note that our corre- 
lation length exponent is close to the v = 1.6 found by 
Onoda et al. [30] in a model of a metal to TI transition. 
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